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PREFACE 

This   investigation  was   oonduoted  at  the  Ifeited  States  Naval  Post- 
graduate Sohool,  Monterey,   California  as  the   thesis  requirement  for   the 
degree   of  Master  of     Science   in  Aerology* 

For  help  and  advice  received  in  its  preparation  the  author   is   indebted 
to  Associate  Professor  G*  J«  Haltiner   of  the  Postgraduate  Sohool. 
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I.      INTRODUCTION 

In  the  attempt  to  subjeot  the   state  of  the  atmosphere  to  a  mathe- 
matioal  analysis,   and   in  particular,   to  explain  the  growth  and  behavior 
of  wave   like  disturbances   in  the  westerlies,   the  meteorologist   is  con- 
fronted with  two  major  problems: 

(a)  An   imperfect  knowledge  of  the   instantaneous  state   of  the 
atmosphere  and 

(b)  The  mathematical  difficulties   inherent  in  the  attempt  to 
integrate  the  equations  of  motion* 

Adequate   observational  data   is   available  to  permit  us  to  draw  some 
conclusions  as  to  the  "average**   conditions   in  the  atmosphere,   or,    in  the 
case   of  choosing  a  particular  model  to  analyse,   a  '^plausible"   condition* 

The  mathematical  problem  is   simplified  by  considering  the  effect  of 
a  small  perturbation  on  a  known  zonal  flow  and  neglecting  all  differentials 
in  the  resulting  eqiiations  of  higher  order  than  the  first*     This  results   in 
a  set  of  linear  differential  equations   —  the  perturbati<Ma  equations*     The 
mathematical  details  have,    in  most  cases,   been  further  reduced  by  making 
certain  simplifying  assvsnptions   in  the  models*     Although  such  assumptions 
introduce  artificialities  when  compared  with  the  atmosphere,   they  essent- 
ially resolve  the  problem  to  the  study  of  dynamic   instability  — •  that 
portion  of  the  stability  or   instability  contributed  by  the  nature   of  the 
flow  pattern   itself*     ICuo     [4 J      states  that,    in  studying  two-dimensional 
nondivergent  flow. 
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Similar  results  would  be  obtained  for  the  very  general  case 
with  vertical  motion,    solenoids,   and  divergence  and  convergence, 
if  certain  approximations  are  made   in  reducing  the  partial   dif- 
ferential equations  for  a   certain  dependent  variable  to  an  ordi- 
nary differential  equation* 

The  technique  of  using  the  perturbation  equations   is  not  new,   having 
been  used  over  fifty  years  ago  by  Helmholti,  Rayleigh  and  others.     In  1928 
Solberg  showed  that  unstable  waves  can  develop  in  a   sloping  surface  of  dis- 
continuity between  two  parallel  streams  of  different  density,    supporting 
the  results  of  his   synoptic   studies  of  wave  development  on  the  polar  front* 
In  1939  Rossby  developed  his  widely  known  equation  for  the  speed  of  long 
waves   in  the  westerlies   in  the  special  case  of  constant  tonal  motion  in  a 
homogeneous  incompressible  atmosphere.     Since  that  time,  a  great  number  of 
meteorologists  have  applied  the  theory  to  a  variety  of  models*     Machta      [s]    , 
using  Rossby* s  model,   assumed  a   solution  for  the  perturbation  stream  func- 
tion which  allowed  for  a  tilt  to  the  trough  line«     Haurwitz  also  used  Rossby* s 
model  but  allowed  for  the  curvature  of  the  earth  and  obtained  only  slight 
modification  to  Rossby* s  results.     Charney     [^ij   considered  a  special  case  of 
baroclinio  three-dimensional  flow« 

Kuo     [4]   developed  general  stability  criteria  for  a  two-dimensional, 
barotropic  nondivergent  atmosphere   using  general,   symmetric  velocity  profiles* 
Eavirwitz      [3j     considered  pure  shearing  waves  between  two  constant  sonal 
two-dimensional  wind  fields  on  a  rotating  planet.     The   case  of  a   shear   zone 
of  finite  width  was  considered  neglecting  the  variation  of  the  ooriolis  para- 
meter with  latitude.     This  problem  in  considered  briefly  in  this   paper. 

It  may  be  noted  that  in  nearly  all  of  the  papers   studied  only  a  portion 
of  the  c(Mnplete  problem  was   investigated.     Some  are   concerned  with  the  question 
of  whether   or  not  a  perturbation  will  be  unstable.     Others  have  eliminated 
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potentially  unstable  waves  from  their  solutions  (as  will  be  done  in  this 
paper)  and  proceed  to  establish  the  speed  with  whioh  the  long  stable  waves 
will  move.  This  approach  is  entirely  reasonable  because  if  the  wave  is  un- 
stable, the  solution  is  valid  only  in  the  initial  development  and  the  sig- 
nifloanoe  of  the  real  values  of  wave  velocity  is  limited  to  this  time.  Also, 
as  Charney  tl]  points  out,  the  evaluation  of  the  complex  roots  may  present 
great  difficulties.  On  the  other  hand,  if  the  wave  is  essentially  stable, 
the  second  order  differentials  whioh  have  been  neglected  in  the  pertur- 
bation equations  will  continue  to  be  small  and  the  value  of  wave  velocity 
obtains  is  valid  for  an  appreoiable  length  of  time* 

In  this  paper,  only  the  long,  slow  moving  or  retrograding  waves  are 
considered*  The  frequency  equation  resulting  from  the  solution  of  the 
boundary  systems  determined  by  a  sinusoidal  velocity  profile  and  a  linear 
shear  tone  are  determined.  The  frequency  equation  for  a  probability  curve 
profile  was  not  fully  developed  due  to  mathematical  difficulties  encountered 
in  the  application  of  one  of  the  boundary  conditions.  This  matter  is  dis- 
cussed in  some  detail. 

The  roots  of  the  frequency  equation  in  the  case  of  the  sinusoidal 
velocity  profile  are  evaluated  (for  a  particular  wave  length)  for  a  synoptic 
situation  in  which  the  horizontal  wind  velocity  profile  was  in  fair  agree- 
ment with  the  model.  The  results  are  found  to  be  of  the  right  order  of 
magnitude  and  show  that  there  may  be  more  than  one  possible  value  of  wave 
speed  under  a  given  set  of  conditions. 
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He  shall  now  consider  three  velocity  profiles  and  derive  the  equations 
relating  wave   speed  and  wave  number.     The   perturbation  equations   to  be   used, 
as  developed  by  Haurwiti      [2l     ,   assume   incompressible  flow.     We   shall  further 
assume  two-dimensional  nondivergent  flow.     The  variation   of  the  coriolis 
parameter  with  latitude   is  assumed  to  be  a   linear   function  of  y.      (Appen- 
dix II).     Under  these  conditions,    the  perturbation  equations  reduce  to 


'"  {t-u)     -h    V  (  U"-  3    --^^  C    -k-  o(^  u)    -    0 


(1) 


as  derived  in  Appendix  I« 

A  complete  analysis  of  this  equation  must  inoliade  a  stvKiy  of  the  be- 
havior of  V  in  the  neighborhood  of  the  singular  point  c  s  U.  Such  a  point 
must  occur  if  the  wave  is  moving  with  a  speed  between  Uq  and  Uj^^j^.  Kuo  [4] 
has  investigated  this  point  for  the  case  of  general  symmetric  velocity  profiles 
by  including  friotional  terms  in  the  above  equation.  These  terms  are  small  com' 
pared  with  the  others  except  when  o  -  U  is  small. 

He  shows  that  if  a  singular  point  exists,  a  stable  wave  can  exist  only 
if  0' -^  =  o     at  some  point  on  the  velocity  profile.   If  this  condition  is 
met,  the  speed  of  such  a  stable  wave  is  the  value  of  U  at  that  point.  Kuo 
calls  this  the  oritical  velocity*  Any  wave  moving  with  less  than  this 
velocity  (but  greater  than  Uq)  will  be  unstable,  while  those  moving  with  a 
velocity  greater  than  the  critical  velocity  will  be  damp«d.  He  further  shows 
that  for  the  long  slow  moving  waves  (c  <,Uq),  the  amplifying  or  dampening 
factors  are  small  or  eero  and  such  waves  are  essentially  neutral. 
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In  the  solutions  to  be  disoussed,  no  closed  form  was  found  and  the 
power  series  solutions  obtained  are  valid  only  for  (o  -  U)  ^  0,  the  slow 
moving  or  retrograding  waves.* 

*  Kuo  ([4I  has  shown  that  the  minimum  wave  length  of  suoh  waves  is  about 
SOO  km  under  typioal  shear  conditions* 
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II.      SINUSOIDAL  VELOCITY  PROFILE 

The  first  problem  to  be  considered  is  the  solution  of  equation   (l) 
where 

U    =     U,     -»■    U,     lU^yy.     X-H  O      ^  H       ^      "^ 


(2) 


A  series  solution  for  v  in  the  sinusoidal  region   is  obtained   in  the  oon* 
rentional  manner  as  sho?m  in  Appendix  III,     It   involres  two  arbitrary 
constants  which  must  be  eliminated  by  the   imposition  of  two  boxzndary 
conditions* 

Euo     [4]    assumes  that  for  symmetric  velocity  profiles  and   symmetric 
boundary  conditions,   there  will  be  either  symmetry  or  antisymnwtry  in  the 
amplitude  of  the  perturbation  stream  fimction,   and  therefore   only  half  the 
sone  need  be   considered*     He  then  points  out  that  it  would  be   unlikely  to 
have  ▼  «  0  along  the   line  of  maximum  wind,   thus  reducing  the   study  to  one 
of  symmetric  disturbances* 

Kuo's  development   is   in  terms   of    ^   ,   the  perturbation  stream  function 
instead  of  v*     It  can  be  shown  that,   under   the  assumptions    involved,    equa- 
tion  (1)   in  either  variable   is   identical* 

The  above   statements  result   in  the   conclusion  that  v  has  a  maximum 
where  U  is  a  maximum;   that   is, 

(In  this  and   subsequent  work,  the  quantity  in  parentheses  refers  to  the 
point  at  which  the   function   is  to  be  evaluated) 
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This  equality  is  used  as  the  first  boundary  condition  and  results  in 
the  elimination  of  one  of  the  arbitrary  constants. 

Where  y  <  0,  in  the  field  of  constant  tonal  wind,  it  can  be  shown 
that  the  solution  of  (1)   is 

V  =  K  e 

Then  , 

V     -  A.  v  =  o 

This   condition  must  be  true  at  all  points   in  the   Eone,   therefore   it   is 
true  at  y  s  0*     Since  y  s  o  is  common  to  both  this  zone  and  the  sinusoidal 
tone,  the  statement   is  true  at  this  point  when  applied  to  the   sinusoidal 
zone*     The  second  boundary  condition  therefore   is 

The   imposition  of  the  boundary  conditions   is  shown  in  Appendix  III 
and  results   in  a  fourth  degree  equation  in  d,  where  d   s  c   -  U   .     The   four 
roots  of  this  equation  may  be  evaluated,   by  rather  tedious  methods,    after 
selection  of  the  physical  parameters   involved* 

Two  synoptic  situations  in  which  the  velocity  profiles  at  the  600  mb. 
level  were  in  fairly  good  agreement  with  that  used  in  this  development  are 
considered  as  examples*  These  are  illustrated  on  the  following  page.  The 
parameters  of  the  upper  profile  were  evaluated  and  the  equation  determined 
for  a  wave   length  of   6000  km.     This  resulted   in 

(C  -  H  S'  Ji     -    1  8  10  4*  -  UE  ^0  4.  -  51  00    ^   o 
The  roots   of  this  equation  are   (approximately) 

d  «  -2  m/sec,    -1.8  nv/sec,    -22  m/seo*,    71*5  n^/seo . 
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The  last  root  may  be  discarded  since  we  are  not  considering  waves  which 
move  faster  than  minimum  wind  velocity*  This  gives  the  possible  values 
of  wave  speed  as  13*4  m/seo*«  13*6  m/seo.,  and  -6*6  cy^seo*  The  wave  is 
stable  since  all  the  roots  are  real* 
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III.   PROBABILITY  CURVE  VELOCITY  PROFILE 

The  second  velocity  field  to  be  considered  in  the  solution  of  (l) 

is 

-ky" 

(3) 

By  varying  the  parameters,  this  profile  may  be  adjusted  to  give  a  good  fit 
to  a  large  number  of  symmetric  profiles  in  the  vicinity  of  the  jet.  The 
differential  equation  is  solved  by  the  series  method  in  Appendix  IV.  The 
first  boundary  condition  is  the  same  as  in  the  previous  problem  and  when 
applied  results  in  a  solution  for  v  involving  one  arbitrary  constant* 
The  second  boundary  condition,  as  before,  is 

V*  -  rv  =  0. 

However,  the  "point"  at  which  this  must  be  applied  is  at  infinity.  The 

use  of  this  condition  therefore  requires  a  transformation  of  the  equation 

in  which  y  r  ^  ,  The  point  at  infinity  is  then  transformed  to  the  origin 
s 

and  the  boundary  condition  may  be  applied.  However,  the  mathematics  here 
is  rather  complicated  and  the  actual  expression  for  c  in  terms  of  o(  was 
not  obtained* 
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17.   LINEAR  SHEAR  ZCNE  VELOCITY  PROFILE 

The  final  problem  is  oonoernod  with  the  solution  of  (1)  for  a  wind 
profile  which  is  defined  by  the  equations 

u  ^  u^     -i.    ^"  ~  "^V-  ^        o  <.  -j  ^  1 

(4) 
0-0,  j  ^  ^ 

This  model  differs  from  the  previous  two  and  from  those  considered  by  Kuo 
in  that  it  is  not  symmetrioal*  It  is  also  obviously  less  realistic.  The 
solution  for  v  in  the  shear  zone  is  given  in  Appendix  V, 

By  similar  reasoning  as   in  the  proceeding  problems,   the  boundary  con- 
ditions are  found  to  be 

v'fo)     -    A.,  V  (o  )      :^      o 

where  r-,  refers  to  the  zone  of  velocity  U^;  and  r^^  "to  "the  rone  of  ^i\» 

By  applying  these  boundary  conditions  as  indicated  in  Appendix  V,  a 
rather  complicated  equation  containing  an  implicit  relationship  between  c 
and  o<^  results*  The  difficulties  of  finding  the  roots  of  this  type  equation 
are  so  great  that  no  attempt  was  made  to  obtain  numerical  results  for  a  par- 
ticular model. 
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V.     C(»iCLUSICKS 

As  with  any  second  order  differential  equation  with  variable   coef- 
ficients,  the  equation   solved   in  this  paper  does  not,    in  general,    yield 
a   solution  in  terms  of  elementary  functions.     It  is  quite  possible  that 
the  problems  considered  here  do  possess  such  solutions  but  none  was  found. 
An  attempt  was  made   in  each  case  to  transform  the  equation  so  as  to  obtain 
a   solution  in  terms  of  Bessel  or  Lagendre  functions,  with  no  success.     An 
arbitrary  selection  of  -wind  profile  will,    in  general,   be  met  by  the  same 
difficulties. 

The  perturbation  equations  are,  therefore,  readily  integrable  only 
if  the  velocity  profile  is  chosen  specifically  to  achieve  this  purpose. 
A  particular  but  arbitrary  profile  such  as  one  of  those  chosen  for  this 
paper  will,  in  general,  yield  an  equation  in  which  the  algebraic  diffi- 
culties in  obtaining  numerical  results  are  of  such  magnitude  as  to  make 
them  practically  worthless.  These  difficulties  would  be  increased  in  the 
case  of  complex  roots. 

The  restrictive  assumptions  made  in  all  published  studies  in  this 
field  have  been  justified,  to  a  certain  extent,  but  must  not  be  forgotten 
when  attempting  to  apply  results  to  the  atmosphere.  The  processess  of  the 
real  atmosphere  are  much  more  complex  than  we  have  assumed.  The  three 
dimensional  nature  of  the  atmosphere,  divergence,  vertical  motion,  and 
n on -homogeneity  certainly  play  an  important  role,  and  often  a  dominant 
one   in  the  development  and  propagation  of  disturbances  in  the  zonal  flow. 


(12) 


e?»cfc   tritrjR'T 


1  i^UJ     Vi,!..;  '  '.f.j 


'.Oii'Ux- 


i^rj.'itt    i.;i    uJiw    .rvirtj  ^'iij.: 


»i'";t8 


rfoje©  jtI    Ei.«rT 


\?   lU    ;v:,  Ji.i' 


V  ^.  .-iu     ->,(>  «ni.''>:rri     "'"Li,  bf 


o«  no 


ftld!*   •?? 


,11/ 1^ 


.  j;ri»aT' 


.?H:t 


-83-OOt  .'O    10 


JO  '^7x©n*goaii>n»coa 


baa 


BIBLIOGRAPHY 


1«       Charney,   J.  G»     The  Dynamios  of  Long  Yliaves   in  a    Barocllaio  Westerly 
Current.  Journal  of  Meteorology.   4:136-162,   October   1947. 

Zm       Haurwitt,    B-     Dynamic  Meteorology.     New  York,  McGraw-Hill,    1941 

3,       Haurwitt,    B»     The   Instability  of  Wind  Discontinuities  and  Shear 
Zones  in  Planetary  Atmospheres*     Journal  of  Meteorology* 
6x200-206,   June   1949. 

4*       Kuo,   H*     Djnaamio   Instability  of  Two-Dimensional  Non-divergent 
Flow  in  a   Barotropio  Atmosphere.     Journal  of  Meteorology. 
6:106-122,   April  1949. 

6.       Machta,   L«     Dynamic  Characteristics  of  a  Tilted-Trough  Model. 
Journal  of  Meteorology.   6:261-266,   August  1949. 


(13) 


lo  x;txXJ:cf*>.  ,ottl 


v3Cic'iooc;-'0M  'ic   i&:ii^o\.      «y- 


.  ic'Co^i  iiSi-". 


(El; 


APFEMDIX  I. 


DERIVATION   OF  SECOND   ORDER  ORDINARY  DIFFERENTIAL  EQUATION 
IN  V  AND  y  FROi  PERTURBATION  EQUATIONS 


The  perturbation  equations  for  two-dimensional  nondivergent  flow 
in  an   incompressible  fluid  where  the  undisturbed  flow  is  tonal  are: 

^*  ''^  «l;j  f    <Jx  (5) 


^\/ 


^  ^-^^=--^4^ 


at  dx  "J-  ^  ' 

J-!±    ^    ±JL      ^  o  (7) 

In  order  to  study  the  velocity  of  wave  perturbations  in  the  x  direction, 
we  may  assimie  the  following  form  of  solution} 


c«  fir-ct) 
t«r  (y-ct) 


(8) 
(9) 
(10) 


Then  from  (7)  and   (8), 

^  =r  jj-  (11) 

Substituting  (11)  and   (lO)   in   (6)  we  obtain 
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Taking  the  partial  derivative  of  (12)  with  respect  to  y  and  equating  it 
to  the  left  hand  member  of  (6),  using  (9)  and  (11),  we  get  an  equality  in 
which  each  term  contains  the  coefficient  M  •  Dividing  through 

by  this  quantity  reduces  the  equation  to  one  in  which  the  only  dependent 
variable  is  the  amplitude  factor  of  v,  v(y).  For  sake  of  brevity,  v(y) 
will  now  be  designated  as  simply  v.  The  differential  equation  then 
reduces  to 
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This  may  be  considered  as  an  ordinary  differential  equation,  but  it  must 
be  remembered  that  the  final  solution  of  v  contains  the  periodicity  factor 
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APPENDIX  II. 
VARUTION   OP  CORIOLIS   PARAMETER  WITH  LATITUDE 
Let  the   ooriolis   parameter  be  represented  by 


where  9     ie  assumed  constant   in  the   zone  imder  oonsideration*       Then       5=   -^ 

but       V  -     2.  '-*^    A*-v»     <t> 
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and       ^     -    -2-    cu   co^  <i>  ^ 


From  the   sketoh  we   see  that 


or 
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Substituting  this  value  of  ~-      in  (A)  we  have 
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APfENDIX   III. 

DEVELOaftENT   OF  FREQUENCY  EQUATION  FOR  SINUSOIDAL  VELOCITY  PROFILE 

We  shall  now  solve  the   equation  obtained  in  Appendix  ^p*  for  the  case 
of  a  sinusoidal  velocity  profile* 

n 

(13) 
The  sinusoidal  profile   is  bounded  on  both  sides  by  a  field  of  constant 
velocity,   Uq» 

Substituting   (13)    in   (l),  we  have  the  following  equation  to  be 
solved: 

v"(  d   u,  oo^  i.^)  ^\/Cd»5  '^  ^/  +  ^)  .0  (14) 

where 

cl    ^     C  -    0  ^ 

Expressing  sin  by  in  series  form 
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(15) 


and  assuming  a  power  series  solution  for  v, 

V  --   ^.   -  a.y    ^  a./  V-  --  --  ■  ^^^^ 

we  may  substitute   (15)  and   (16)   in   (14)  and  obtain  an  equation  in  positive 

powers   of  y  equal  to  rero.     Upon  equating  the  coefficients  of  each  power 

of  y  to  sero,  we  find  that 
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This  expression  for  ▼  involves  two  arbitrary  constants  which  may  be  elim- 
inated by  the  imposition  of  two  boundary  conditions*  The  first  condition 
to  be  applied,  as  shown  on  page  6,  is  that  v  have  a  maximum  where  U  is  a 
a  maximum,  or. 

Taking  the  first  derivative  of  (17)  with  respect  to  y,  setting   y 
and  equating  the  result  to  lero,  we  have 
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and  upon  solving  for  a^  we  see  that 
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The  resulting  solution  for  v,  neglecting,  for  the  moment,  the  third  and 
higher  powers  of  y,  is 
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By  imposing  the  second  boundary  condition  (see  page   7), 
v(  a)  -     A,   1/  (  0  )      -    0 


we  obtain 
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This  step  shows  that  only  the  first  two  terms  of  the  series  solution  for  v 
are  significant  in  the  development  of  the  frequency  eqxjation  so  the  powers 
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of  y  higher  than  the  first  could  have  been  ignored  after  having  deter- 
mined equation  (I?)*  However,  it  should  be  noted  that  the  solution  of 
a,  in  equation  (IJ)   is  approximate  since  only  the  first  four  terms  of 
the  general  solution  for  v  were  used  in  its  determination*   It  follows 
that  the  frequency  equation  is  approximate. 

Squaring  equation  (21)  and  collecting  the  coefficients  of  powers 
of  d  we  have 
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APPENDIX  IV. 


DEVELOBJENT  OF  FREQUENCY  EQUATION  FOR  PROBABILITY 
CURVE  VELOCITY  PROFILE 


In  this  section,  the  differential  equation 

is  to  be  solved  for  the  oase  in  which  the  undisturbed  wind  velocity 
profile  is  defined  by  the  probability  curve 

0  =  O^  -V  O,  A  (23) 

A  power  series  solution  for  v  is  assumed, 

v=  ag  ^A-.y  -"^^y  -^ (24) 

and  X     is  expressed  in  series  form. 


(2  5) 


The  series  expressions  for  v  and  U  are  then  substituted  in  (1).  The 
coefficients  of  each  power  of  y  are  equated  to  zero  and  we  may  now  obtain 
expressions  for  any  a^,  in  terms  of  a  and  aii 
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a,- 4, where  V^- O.-t-U,     and 

a  and  a.  are  arbitrary  constants.  The  general  solution  for  v  is  therefore 
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We  may  eliminate  one  of  the  constants  by  the  imposition  of  the  boundary 
condition 

The  use  of  this  condition  for  symmetrical  velocity  profiles  is  discussed 
on  page  6. 

Taking  the  first  derivative  of  (26)  with  respect  to  y,  and  setting  y  *  0, 
we  obtain 

a,  =  o 

Equation  (26)  then  reduces  to 


v^  a 


As  mentioned  on  page  10 «  the  application  of  the  second  boundary 
condition, 

leads  to  mathematical  difficulties  and  the  elimination  of  a  and  consequent 
development  of  the  frequency  equation  was  not  accomplished* 
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APPENDIX  V. 
DBVELOHilENT  OF  FREQUENCY  EQUATION  FOR  A  LINEAR  SHEAR  ZCRE 
In  this  section  the  differential  equation 


(1) 


is  solved  for  the  oase  in  which  the  undisturbed  irind  velocity  profile  is 
a  linear  shear  zone.  The  arbitrary  constants  are  eliminated  and  the  fre- 
quency equation  developed  for  the  case  in  which  the  shear  tone  is  termi- 
nated, on  each  side  by  a  field  of  constant  sonal  wind. 
The  shear  tone  is  defined  by 

^  (27) 

A  power  series  solution  for  v  is  assumed 

and  this  series  is  substituted  in  (1)*  The  coefficients  of  each  power 
of  y  are  equated  to  tero  and  the  following  relationships  obtained x 
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As  explained  on  page  11,  the  boundary  conditions  to  be  imposed  are 

\/'  I  u)   .    A  ,  V  ( 0)      =  0  (29) 
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By  substituting  (28)  in  (29)  and  letting  y  «  0  we  find  that 

a,  -  ^  .  a  0 

Equation    (28)  then  reduces   to 
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To  impose  the  second  boundary  condition,  we  substitute  (31)  in  (30)  and 
set  y  = -^  .  After  dividing  through  by  the  common  factor  a^,  we  obtain  the 
frequency  equation  which  becomes,  upon  rearranging: 
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This  is  approximate  since  fourth  and  higher  powers  of  y  have  been  neglected < 
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